Hematopoietic-specific heterozygous loss of Dnmt3a exacerbates colitis-associated colon cancer

A preclinical mouse model demonstrates that genetic alterations in the blood system characteristic of clonal hematopoiesis (CH) contribute to an aggressive solid tumor phenotype. It further identifies cancer angiogenesis as a potential therapeutic target to mitigate adverse CH effects.

Colon cancer is one of the leading causes of cancer-related deaths in developed countries (Ahmed, 2020;Keum and Giovannucci, 2019;Siegel et al., 2019). Inflammatory bowel disease (IBD), including UC and Crohn's disease, is a wellknown risk factor for colon cancer (Beaugerie and Itzkowitz, 2015). Colitis-associated colon cancer (CAC) accounts for 15% of overall mortality among all IBD patients. Compared with sporadic colorectal cancer, patients with CAC often present with multifocal tumors arising from precancerous lesions that are challenging to detect and remove endoscopically, and tend to rapidly develop chemoresistance (Yaeger et al., 2020). Over 20% of patients with colon cancer in the publicly available Memorial Sloan Kettering Cancer Center (MSKCC) clinical sequencing database have detectable CH (Cerami et al., 2012;Gao et al., 2013), which is notably more prevalent than in agematched cancer-free population. Despite screening and lifestyle interventions, most patients present with advanced disease associated with poor outcome. To inform the choice of optimal therapeutic approaches and improve survival rates, better understanding of disease-modifying factors that contribute to the aggressive tumor phenotype is critically needed.
Potential involvement of CH in the pathogenesis of coincident solid tumors has a far-reaching translational impact, yet a better understanding of this relationship is hindered by the lack of animal models. Here, we combined a well-established induction of CAC (De Robertis et al., 2011;McIntyre et al., 2015;Parang et al., 2016) with a BM transplantation (BMT) approach to experimental CH driven by heterozygous Dnmt3a loss (Guryanova et al., 2016a(Guryanova et al., , 2016b. This unique tool enables rigorous interrogation into the role of CH in the pathogenesis of solid tumors driven by unrelated genetic alterations, uncoupled from environmental, lifestyle, iatrogenic, or other confounding factors (Bolton et al., 2020;Coombs et al., 2017;Dawoud et al., 2020;Gillis et al., 2017;Hong et al., 2022;Hormaechea-Agulla et al., 2021;Hsu et al., 2018;King et al., 2020). At the same time, genetic alterations in BM-derived cells are likely to affect their function and contribute to disease conditions outside of the hematopoietic system. While attempts to unveil the role of CH in non-malignant conditions such as atherosclerosis or gout were made (Agrawal et al., 2022;Fuster et al., 2017;Jaiswal et al., 2017;Rauch et al., 2018;Sano et al., 2018), similar studies in cancer models are limited. In this proof-of-concept study, we report that heterozygous inactivation of Dnmt3a restricted to the hematopoietic compartment exacerbates the development of colon cancer phenotype in an inflammation-induced mouse model. We further credential an anti-angiogenesis drug axitinib as a potential therapeutic mitigation strategy.

Results and discussion
Heterozygous loss of Dnmt3a in the BM leads to accentuated CAC phenotype To explore the relationship between genetic alterations in the blood system and the severity of coincident CAC, we examined presence of CH mutations in paired blood and tumor samples in a cohort of 66 patients with CAC treated at MSKCC. We observed a strong trend toward a more advanced disease at diagnosis among patients with detectable CH by MSKCC-IMPACT (Cerami et al., 2012;Gao et al., 2013;Fig. 1 A, P = 0.056, one-sided Fisher's exact test comparing patients with local [stages 0, 1, and 2] and advanced [stages 3 and 4] disease, and Table 1). Patient age at diagnosis could not explain the difference in CH prevalence ( Fig. 1 B), consistent with previous reports implicating that the impact of CH on solid tumor outcomes was not due to differences in age. There was no correlation between the size of CH clone measured as the variant allele frequency and CAC stage (P = 0.51, Mann-Whitney test) or prior chemotherapy exposure (P = 0.27, Mann-Whitney test), bearing in mind that the cohort was too small to lend sufficient statistical power for such analyses, and CH mutations were only seen in patients with advanced stage disease (3 and 4). To directly test the causal relationship between presence of CH and CAC pathogenesis, we performed a proof-of-concept animal study.
To experimentally model CH in mice, we adopted a BMTbased approach using Dnmt3a +/f :Mx1-Cre + and Dnmt3a +/+ : Mx1-Cre + (WT control) mice as donors after successful Crerecombination induced by intraperitoneal (IP) administration of poly(I:C), yielding Dnmt3a +/Δ and Dnmt3a +/+ mice. Fully engrafted Dnmt3a +/Δ and control Dnmt3a +/+ BM chimeras (Fig. S1 A) underwent a CAC induction protocol initiated by a single injection of azoxymethane (AOM, 10 mg/kg, IP) followed by three cycles of dextran sulfate sodium (DSS) salt exposure (3% in sterile drinking water; Santiago et al., 2019;Tanaka et al., 2003;Fig. 1, C and D). AOM/DSS CAC is a gold-standard autochthonous model with high penetrance and predictable latency (De Robertis et al., 2011;McIntyre et al., 2015;Parang et al., 2016). The use of immunocompetent animals preserves all modes of tumor-microenvironment interactions. AOM/DSS chemically induced CAC produces a spectrum of mutations that captures genetic landscapes observed in human colon cancer (De Robertis et al., 2011) and recapitulates human pathology. This unique tool enables rigorous interrogation into CH involvement in the pathogenesis of solid tumors driven by unrelated genetic alterations, which can be readily extended to other cancer types. Colonoscopy at 9 wk after initiation of the AOM/DSS treatment demonstrated increased colon wall opacity, visible bleeding, and numerous fibrin patches, indicating heightened colon pathology, and detected more, larger tumors in accessible colon regions in Dnmt3a +/Δ BM chimeras compared with Dnmt3a +/+ -engrafted controls ( Fig. 1 E). Modified murine endoscopic index of CAC severity (MEICS; Fung and Putoczki, 2018) independently scored by two blinded investigators was significantly elevated in the Dnmt3a +/Δ group (Fig. 1 F).
To further investigate the effect of Dnmt3a heterozygous loss in the hematopoietic system characteristic of CH on CAC induction, mice were euthanized for comprehensive colon evaluation 10 wk after initiation of the AOM and DSS treatment. Colon gross pathology ( Fig. 1 G) demonstrated significantly elevated penetrance of the tumor phenotype (incidence of large tumors >2 mm in diameter, Fig. 1 H), higher large tumor number ( Fig. 1 I), and increased total tumor burden per colon ( Fig. 1 J) in Dnmt3a +/Δ chimeras. Together, these findings suggest that hematopoietic-specific Dnmt3a haploinsufficiency promotes both cancer initiation and progression in the context of CAC. To more closely model human CH where mutant clones contribute a fraction of mature cells in the peripheral blood, we next transplanted WT recipient animals with 10% Dnmt3a +/Δ (or CAC patients with detectable CH at diagnosis (n = 9) tend to be more likely to present with an advanced stage of the disease (stage 3 or 4) than those without CH (n = 57; Mann-Whitney test, P = 0.074). (B) Age (yr) of CAC patients with or without CH at diagnosis (n = 57 or 9, Mann-Whitney test, P = 0.9). (C and D) Experimental workflow using a mouse model of BMT-based CH and AOM/DSS-induced CAC (C) and timeline (D). (E and F) Representative colonoscopy findings by the Tele Pack Vet X miniendoscopic system in a mouse model of CH (E) and corresponding MEICS scores in Dnmt3a +/+ and Dnmt3a +/Δ BM chimeras (F, n = 14-13 from two replicate experiments, Student's t test, P = 0.05). The diameter of the endoscopy probe is 1.9 mm. (G-J) Representative gross pathology of colons with tumors (black arrowheads, G, bar, 10 mm), proportion of animals presenting with tumors (>2 mm in diameter) in nine replicate experiments (H, Mann-Whitney test, P = 0.019), number of large tumors (>2 mm in diameter) per animal (I, n = 55-53, Mann-Whitney test, P = 0.0015), and cumulative tumor burden per animal (J, n = 55-53, Mann-Whitney test, P = 0.0017). (K) Mice transplanted with a mixture of 10% Dnmt3a +/Δ and 90% WT competitor BM (10% CH BMT) develop more tumors following CAC induction by AOM/DSS, compared with 10% Dnmt3a +/+ controls (n = 14/group in one experiment, Mann-Whitney test, P = 0.058). *, P ≤ 0.05; **, P ≤ 0.01; ns, not significant.     Dnmt3a +/+ control) mixed with 90% WT competitor BM cells distinguishable by CD45.2 and CD45.1 pan leukocytic markers. Following the AOM/DSS CAC induction protocol, we observed a higher number of colon tumors in the Dnmt3a +/Δ -CH 10% BM chimeras ( Fig. 1 K) and a higher tumor tropism of CH-derived myeloid cells (Fig. S1 B), demonstrating that a small Dnmt3a-CH clone was sufficient to promote cancer initiation.
Dnmt3a haploinsufficiency in the BM leads to accentuated pathological features of CAC Given that full Dnmt3a +/Δ BM chimeras have higher tumor burden and larger tumor size compared to WT controls, we performed histopathology analysis on colons from both groups based on a modified quantitative scoring system for DSSinduced murine CAC (Arthur et al., 2012;Cooper et al., 1993;Dieleman et al., 1998;Karrasch et al., 2007). H&E-stained paraffin sections of Swiss-rolled colons showed marked immune infiltration, extensive ulceration and dysplasia of colonic epithelium, and more frequent adenocarcinoma formation with occasional submucosal invasion in Dnmt3a +/Δ -reconstituted animals ( Fig. 2 A). In comparison, WT control chimeras showed moderate dysplasia, fewer adenoma polyps, mild ulceration, and more tissue regeneration. Histology scores based on four parameters (immune infiltration [0-3], ulceration [0-3], morphology of colonic epithelium [0-4], and neoplasms [0-4]), independently assigned by two blinded investigators, were significantly higher in the Dnmt3a +/Δ group than in WT (Fig. 2, B and C). We found that in Dnmt3a +/Δ BM chimeras, increased proliferation of colonic epithelium occurred early in the carcinogenesis process, evidenced by a higher proportion of Ki67marked cells per crypt in the mucosal layer after the first cycle of DSS (Fig. 2, D and E). Overall, Dnmt3a haploinsufficiency in the BM yields more advanced CAC histopathology, consistent with a more severe tumor burden.
Gene expression profiling identifies signatures of accentuated CAC tumorigenesis in animals with experimental CH driven by Dnmt3a +/Δ To identify specific molecular mechanisms likely driving accentuated CAC tumorigenesis in animals with Dnmt3a +/Δ hematopoiesis, we profiled tumor transcriptomes from WT and Dnmt3a +/Δ chimeras by bulk RNA sequencing (RNA-seq; Table  S1). We identified 297 differentially expressed genes (log 2 fold change > 1, P-adjusted < 0.05, Fig. 3 A) with most being upregulated rather than downregulated in the Dnmt3a +/Δ group (241 up and 56 down, Table S2). Gene set enrichment analysis (GSEA; Subramanian et al., 2005) using the HALLMARK collection of gene signatures (Liberzon et al., 2015; Fig. 3 B) detected significant positive enrichment of cancer-and proliferation-related pathways such as epithelial-mesenchymal transition, Wnt/ β-catenin, and angiogenesis (Fig. 3, C-E), along with E2F, mitotic spindle, G2/M checkpoint, and MYC signaling, and shifts in metabolism such as negative enrichment of oxidative phosphorylation (Fig. 3 F) and adipogenesis/fatty acid metabolism (Table  S3). Elevated expression of the Wnt/β-catenin pathway genes is notable (Fig. 3 C) given its known role in driving colon cancer and association with poor prognosis (Jackstadt et al., 2020), which most likely reflects a more advanced stage of the disease evidenced by increased tumor size and overall cancer burden in Dnmt3a +/Δ -transplanted animals (Fig. 1, G and J). The leading edge of the angiogenesis signature was dominated by vascular endothelial-specific genes along with fibroblast growth factor receptor 1 (Fgfr1), angiogenesis-promoting guanine nucleotide exchange factor for the Rho family of Ras-related GTPases Vav2, and a homeobox transcriptional repressor Msx1 (Table S3), while Vegfa/b or Vegfr1/2/3 were not differentially expressed. Consistently, intratumoral levels of vascular endothelial growth factor (Vegf) did not differ between Dnmt3a +/Δ and control Dnmt3a +/+ BM chimera (Fig. S2 A). Overall, these results are consistent with our findings of increased proliferation within the epithelial layer of regenerating colons in Dnmt3a +/Δ chimeric mice likely leading to a more advanced CAC pathology.
Treatment with an angiogenesis inhibitor eliminates the tumor-promoting effect of Dnmt3a +/Δ BM To further investigate our finding of an enriched angiogenesisrelated gene expression signature detected by RNA-seq, we performed immunofluorescent staining of colon tumors using anti-CD31 antibody that marks endothelia. CAC tumors from animals grafted with Dnmt3a +/Δ BM had increased CD31 staining area compared with Dnmt3a +/+ -reconstituted controls, consistent with more dense vascularization (Fig. 4, A and B) known to promote tumor growth. Hence, we hypothesized that cancer angiogenesis in animals with experimental CH may be targetable therapeutically to mitigate accentuated tumor phenotype. Axitinib is an orally bioavailable small-molecule angiogenesis inhibitor that is FDA-approved for treating advanced renal carcinoma with clinical trials for other cancers ongoing (Gross-Goupil et al., 2018). Axitinib is a tyrosine kinase inhibitor of Vegf receptors −1, −2, and −3 that is also active against plateletderived growth factor receptor (PDGFRβ) and stem cell factor receptor c-KIT/CD117 as part of its potent anti-angiogenic mechanism of action (Escudier and Gore, 2011). After confirming that axitinib potently inhibited cancer angiogenesis in a CAC model (Fig. S2, B and C), we tested if it also abolished the colon tumor-promoting effect of Dnmt3a +/Δ BM (Fig. 4 C). A 10-wk course of axitinib treatment eliminated the difference in vascular density between genotypes (Fig. 4, D and E). At the same time, a substantially lower proportion of Dnmt3a +/Δ -BM grafted animals had large tumors compared with vehicle controls (Fig. 4 F). Treated animals in the Dnmt3a +/Δ -BM group developed dramatically fewer large tumors and presented with significantly lower tumor burden per colon (Fig. 4, G and H), while the tumor phenotype in the Dnmt3a +/+ -grafted mice was only slightly affected by axitinib administration. Since axitinib is a multikinase inhibitor that also binds to the stem cell factor receptor cKit essential to the function of hematopoietic stem and progenitor cells, we wondered if it affected hematopoiesis. While there were no differences in the major lineage composition in the peripheral blood of vehicle-treated mice transplanted with Dnmt3a +/Δ and control Dnmt3a +/+ kept on DSS, administration of axitinib resulted in a decrease in circulating myeloid cells and a coincident increase in T cells in Dnmt3a +/Δ BM chimeras (Fig. S3). At the same time, axitinib treatment reversed expansion of the Lin − Sca1 + cKit + CD48 − CD150 + long-term (LT) HSCs in the BM of Dnmt3a +/Δ -transplanted animals (Fig. 4 I). These findings indicate that in addition to inhibiting tumor angiogenesis, axitinib may also act indirectly through rebalancing dysfunctional Dnmt3a-CH hematopoiesis in the context of inflammation. Overall, axitinib treatment reversed the heightened CAC tumor phenotype in Dnmt3a +/Δ -chimera, suggesting its potential translational utility for mitigating unfavorable effects of CH in cancer patients.
Although in our study using a CAC model Dnmt3a-CH promoted both tumor initiation and progression, consistent with prior pan-cancer retrospective clinical observations (Coombs et al., 2017), this effect is likely context specific. Thus, in patients with metastatic colorectal cancer enrolled in the FIRE-3 clinical trial comparing standard chemotherapy plus cetuximab or bevacizumab, both of which inhibit angiogenesis, presence of CH, and CH-DNMT3A specifically, was associated with extended overall survival (Arends et al., 2022). Consistently, in our clinical cohort, patients with CH presented with more advanced disease at cancer diagnosis, although the study was not powered to evaluate survival. These contrasting observations may reflect disparate effects of DNMT3A mutations in different branches of the immune system within the tumor microenvironment tipping the balance between anti-tumorigenic and pro-resolving phenotypes (Leoni et al., 2017;Prinzing et al., 2021). Similarly, it is possible that different DNMT3A mutation types may exert disparate effects on antitumor immunity. Thus, in our cohort, a patient with DNMT3A(R882) CH did not develop distant metastases during the observation period while both patients with non-R882 mutations did. This observation, while anecdotal, lends credence to a possibility of different CH alterations imparting divergent functional consequences. While outside the scope of this study, these specific mechanisms remain to be investigated in the future. Further, the effects of Dnmt3a-CH are likely pleiotropic and extend beyond anticancer immunity, illustrated by the identification of changes in tumor metabolism and the striking finding of enhanced tumor angiogenesis in animals with experimental CH. In support of this, targeted inhibition of angiogenesis with an FDA-approved small molecule tyrosine kinase inhibitor axitinib abrogated the CAC tumorpromoting effect of experimental Dnmt3a-CH, directly by decreased tumor vascular density and indirectly by normalizing aberrant Dnmt3a-CH hematopoiesis. These results identify an actionable therapeutic strategy to mitigate the tumor-promoting effect of coincident DNMT3A-CH that can be tested in the clinical setting.
Together, these results highlight a significant difference of molecular pathophysiology effected by genetic alterations within the hematopoietic compartment, despite identical mode of CAC induction. It indicates that alteration of Dnmt3a in the BM has a profound impact on molecular pathogenesis of CAC through multiple non-tumor-cell-autonomous mechanisms, some of which may be targetable. These findings, for the first time, solidify the causal relationship between CH and the severity of solid tumors and identify potential therapeutic strategies. Further study is needed to investigate underlying molecular mechanisms including immune involvement.

Patients
Patients with CACs were identified from a database of genomically annotated CAC cases maintained under MSK Institutional Review Board protocols 15-297 and WA0143-14. All participating patients signed informed written consent for matched tumor and normal sequencing (MSK Institutional Review Board protocol 12-245), and next-generation sequencing was performed with the MSK-IMPACT assay (Cheng et al., 2015). MSK-IMPACT is a hybridization capture-based next-generation assay encompassing all exons of >340 genes. It is validated and approved for clinical use by the New York State Department of Health Clinical Laboratory Evaluation Program. The sequencing test utilizes genomic DNA extracted from formalin-fixed paraffin-embedded tumor tissue as well as matched patient blood samples. DNA is sheared and DNA fragments are captured using custom probes. MSK-IMPACT contains most of the commonly reported CH genes with the exception that earlier versions of the panel did not contain PPM1D or SRSF2.

Mice
Animals were housed at the University of Florida Cancer & Genetics Research Complex specific pathogen-free animal facility; all animal studies were approved by the University of Florida Institutional Animal Care and Use Committee under protocol #201909474. A conditional Dnmt3a knock-out line was previously described (Guryanova et al., 2016a;Nguyen et al., 2007). Dnmt3a +/fl and Dnmt3a +/+ mice with Mx1-Cre deleter on a C57BL6/J background were generated via in-house breeding. To achieve inducible hematopoietic-specific excision, 8-12-wk-old Dnmt3a +/fl :Mx1-Cre and littermate WT control mice received five IP injections of poly(I:C) (#tlrl-pic-5; InvivoGen) every other day. Mx1-Cre-driven recombination was validated by PCR using genomic DNA from peripheral blood mononuclear cells 2 wk after the last injection, at which point the BM was harvested and  The Jackson Laboratory) through tail vein injection. Successful engraftment was confirmed by CD45.1/CD45.2 peripheral blood chimerism 2 mo after BMT. Mice of different genotypes were cohoused to control for possible cage effects. Animals received unique IDs for the purpose of blinding; investigators were unblinded after analyses were complete. Peripheral blood was collected by submandibular puncture. Complete blood counts were obtained using HESKA HT5 automated veterinary hematology analyzer. Both male and female recipients were used with equivalent results.

CAC induction
CAC is induced in a 10-wk-long protocol. Successfully engrafted mice received a single dose of AOM (Cat# A4586, 10 mg/kg; Sigma-Aldrich) by IP injection. 1 wk later, DSS (AAJ6360622, MW~40,000; Alfa Aesar) was provided in sterile drinking water at 3% (wt/vol) as the sole source of water during each of the three DSS treatment cycles with 2 wk of sterile drinking water between cycles. As the effect of DSS may vary between mouse strains, genders, and housing facilities, each lot was tested to determine the optimal DSS dosing. Body weights were monitored daily during DSS administration, and animals losing more than 25% body weight were euthanized as past humane endpoints approved by the Institutional Animal Care and Use Committee. Colon microbiota was normalized by mixing cage bedding of all experimental groups every other week to eliminate cage-to-cage variability. Since AOM is a systemic carcinogen, all animals were bled monthly to monitor for hematologic malignancies. Animals with hematopoietic disorders diagnosed by complete blood counts and flow cytometry analysis of peripheral blood mature lineages (Mayle et al., 2013) were excluded as potential confounders. At day 70, animals were bled and euthanized for comprehensive analysis. Flow cytometry data were collected on a 5-laser 16-parameter BD LSR Fortessa instrument and analyzed by FlowJo v10 using the antibody cocktail in Table 2. For in vivo axitinib treatment studies, axitinib (Cat. # S1005; SelleckChem) was diluted in 0.5% carboxymethylcellulose sodium salt to achieve 2.5 mg/ml suspension. Fully engrafted mice were administered 25 mg/kg axitinib three times/wk (Monday-Wednesday-Friday) by oral gavage throughout the duration of the CAC induction. For the analysis of axitinib's impact on hematopoiesis, AOM was omitted to avoid its carcinogenic effect as a potential confounder.

Colonoscopy in live animals
Colonoscopy was performed after second and third DSS cycles using the Tele Pack Vet X mini-endoscopic system (KARL STORZ Veterinary Endoscopy) equipped with a rigid 64301 AA Hopkins Straight Forward Telescope 0°with a 10-cm fiber optic light transmission as previously described (Becker et al., 2006;Lippert et al., 2009;Mottawea et al., 2016;Uronis et al., 2009). Briefly, mice were anesthetized with isoflurane and placed ventral side up on a heating pad. An endoscope was carefully inserted into the rectum up to 4 cm under visual control with slow air flow to keep the colon inflated, and withdrawn slowly while recording the localization and size of colon abnormalities. Colonoscopy videos were independently scored by two blinded investigators using MEICS (Fung and Putoczki, 2018) that combines five scoring criteria (colonic wall thickening [0-3], vascular pattern [0-3], fibrin formation [0-3], stool consistency [0-3], and tumor diameter [0-5], for a total score ranging 0-17).

Tissue dissection and processing
Colons were dissected immediately after euthanasia. After washing out the feces with 2-3 ml sterile PBS, colons were laid on wet Whatman blotting paper and cut longitudinally. Tumors were counted and dimensions in millimeters were measured with digital calipers used to calculate tumor volume V = (Width 2 × Length)/2. For histological analysis, each dissected colon was Swiss-rolled from the distal end.

Feng et al. Journal of Experimental Medicine
Scientific). Frozen samples were cut into 12-µm-thick sections at −20°C, bound to positively charged microscope slides (Cat# SM2575; ASI) at room temperature (RT), and stored at −20°C.
For immunofluorescent staining, slides were thawed at RT for 20-30 min and tissue permeabilized with 0.2% Triton-X 100 in PBS at RT for 10 min, and then thoroughly washed with PBS. After blocking with 10% donkey serum in PBS for 1 h, slides were incubated with goat anti-mouse CD31 primary antibodies (#AF3628; 1:100; R&D Systems) in blocking buffer overnight at 4°C and stained with secondary AlexaFluor633 donkey anti-goat antibodies (#A21082; 1:200; Invitrogen) for 1 h at RT, three washes with PBS after each incubation, and mounted with Prolong Gold Antifade with 49,6-diamidino-2-phenylindole (DAPI; P36935; Invitrogen) as a counterstain for nuclei. Images were acquired with a Leica DMi8 confocal microscope equipped with a DFC7000 camera using a 20× HC PL Fluotar objective (Leica) using uniform settings across all slides. The LAS Navigator function was used to generate a merged image of the whole cross-section. All images were processed identically using Fiji/ ImageJ (National Institutes of Health).
Tumor vascular density was quantified as CD31 (endothelial marker) staining area per mm 2 of each tumor using ImageJ v1.53k. First, tumors were manually selected using the DAPI channel of the multiple-layer .lif files of the colon Swiss roll. Next, each tumor was cropped and tumor area (in mm 2 ) was measured. To reduce background, uniform adjustments were made to the brightness, contrast, and threshold settings in the CD31 channel for all samples. The area of CD31 staining (in mm 2 ) was measured and tumor vascularization was calculated as the ratio of CD31 area to tumor area.

Measurement of intratumoral Vegf concentrations
Tumors excised from surrounding colon tissue were homogenized in 500 μl of 1% NP-40 lysis buffer (150 mM NaCl, 20 mM Tris HCl, pH 7.5) supplemented with protease and phosphatase inhibitors using 0.5 ml of 800 µm glass homogenization beads in screw-top blast microtubes placed in the Benchmark Bead Blaster 24 Microtube Homogenizer at 6 m/s for 30 s with 2-min breaks between the three cycles. Protein concentration in supernatants was quantified using BCA Protein Assay (Life Technologies). All samples were diluted to a uniform protein concentration of 4 mg/ml and submitted to Eve Technologies Corp. for Murine 31-plex cytokine and chemokine discovery arrays.
Analysis of the stem/progenitor compartment in the BM Fully engrafted animals administered with axitinib (or vehicle control) and three cycles of DSS were harvested on day 70. Peripheral blood was collected by a submandibular puncture. BM was isolated from long bones by spin-flush method. Red blood cells were depleted by incubating with ammonium-chloridepotassium lysis buffer (A1049201; Thermo Fisher Scientific) on ice for 15 min. Analysis of the donor-derived (CD45.2 + ) stem and progenitor cell populations was performed on single-cell suspensions by flow cytometry. LT-HSCs were defined as Lineage − Sca1 + c-Kit + (LSK) CD150 + CD48 − , granulocyte/macrophage progenitors-LSK CD16/32 + CD34 + , using the antibody cocktail in Table 3.
Sample preparation for RNA extraction, sequencing, and analysis After dissection, each colon tumor was submerged in 200 μl RNAlater (AM7021; Thermo Fisher Scientific) in 1.5-ml Eppendorf tube, flash-frozen in liquid nitrogen, and stored at −80°C. Tumor tissue was thawed on ice and transferred into a blaster tube containing 500 μl buffer RLTplus (1053393; Qiagen) and 200 μl glass beads (12621-148; VWR), then homogenized using BeadBlaster 24 Microtube Homogenizer (Benchmark) at 6 m/s for 30 s three times, with 2-min breaks for cooling. Total cellular RNA was isolated with RNeasy Microprep kit (74004; QIAGEN) and quality controlled on a 4150 TapeStation (Agilent). RNA samples with an RNA integrity number 7.0 or higher were picked for library preparation and sequencing. RNA was subjected to standard Illumina-based RNAseq library preparation and sequenced on an Illumina NovaSeq6000 using a paired-end 100 bp chemistry to an average depth of 49 million reads/sample. Input sequences were trimmed with Trimmomatic; quality control was performed before and after trimming using FastQC. Retained reads (>96%) were aligned to mouse reference transcriptome mm10/build 38 using STAR 2.7.3a; gene and transcript expressions (raw fragments per kilobase of transcript per million mapped fragments) were quantified using RSEM v1.2.31. From a transcript count matrix, differential gene expression was evaluated with DESeq2 using a log2 fold change cutoff of ±1 and a false discovery rate of 5% (Love et al., 2014). GSEA was conducted using the fgsea package in R. The volcano plot (Fig. 3 A) and heatmap (Fig. 3 C) were produced using the EnhancedVolcano and pheatmap R packages, respectively. All R code used to analyze RNA-seq data will be made publicly available at https:// github.com/RobinsonTroy/CH_CAC_RNAseq.

Statistical analysis
Sample size calculation was based on the primary endpoint: tumor burden (mm 3 ). The number of animals was chosen to ensure 90% power with 5% alpha to detect a difference between groups of one standard deviation (SD) or larger based on variability and technical drop-out rate observed in pilot experiments. All grouped data are presented as mean ± SD. Statistical significance was determined by unpaired parametric Student's t test and by non-parametric Mann-Whitney test after testing for normal distribution. For samples with significantly different variances, Welch's correction was applied. Statistical analyses and visualization of the data were performed using Prism 9.0.2 (GraphPad Inc.). For pair-wise comparisons, P values ≤0.05 were considered significant.

Rigor and reproducibility
For the entire study, both female and male animals were used to control for gender-specific effects. All experimental animals from different genotype groups were co-housed to mitigate potential cage effects. Other variables were kept consistent in all cages and experiments, such as microbiota and DSS quality. Animals developing hematologic malignancies were excluded from analysis as potential confounders. Blinding was achieved by assigning random codes to each animal and sample prior to analysis; investigators were unblinded after results had been recorded.

Data availability
All data underlying Figs. 1, 2, and 4 are available in the published article and its online supplemental material. The data underlying Fig. 3 and Tables S1, S2, and S3 (RNA-seq) are openly available at the NCBI's Gene Expression Omnibus under the accession number GSE213178.  (B) Dnmt3a-CH-derived myeloid cells (CD11b + ) exhibit enhanced tumor tropism in the 10% CH BMT CAC model, as seen by increased CD45.2 chimerism (n = 12-13, one experiment, two-tailed t test with Welch's correction, P = 0.009). **, P ≤ 0.01; ns, not significant. Figure S2. Axitinib suppresses tumor vascularization. (A) Intratumoral concentration of Vegf did not differ between mice transplanted with Dnmt3a +/Δ and control Dnmt3a +/+ BM (n = 14-17 from two independent experiments, Student's t test with Welch's correction, P = 0.55). Vegf concentrations were determined in tumor protein lysates by multiplex ELISA performed by Eve Technology Corp. (B and C) Administration of a small molecule VEGFR and multikinase inhibitor axitinib potently decreases tumor angiogenesis. Mice transplanted with Dnmt3a +/+ BM were treated with axitinib (25 mg/kg per os, three weekly doses) throughout CAC induction. Representative images of colon tumors (B, vasculature detected by endothelial marker CD31 in red; DAPI was used to visualize nuclei in blue; bar, 100 μm). Vascular density (as CD31 staining area per mm 2 ) was quantified at endpoint in 16-18 colon tumors in vehicle and treatment groups from one experiment (C, Student's t test with Welch's correction, P = 0.009). The images of colon tumors from Dnmt3a +/+ chimeras treated with axitinib are also presented in Fig. 4 D. **, P ≤ 0.01; ns, not significant.